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^ ■ Abstract 



Q I The statistical analysis and the spherical wavelet analysis of the SDSS DR7 

■ quasars distribution and of the WMAP CMB anisotropy are performed. We re- 

vealed the qualitative agreement between the angular power spectrum of CMB and 
the angular power spectrum of the quasar distribution on the celestial sphere. The 
angular correlation function and the angular power spectrum of the quasar distribu- 
tion may be described by the power laws oj {'&) ~ i!}^^"^ and ui ~ P^~^ respectively 
with the averaged parameters (7) = 2.08 and (7^) = 1.92. The large quasar groups 
are discovered and they form the fractal set: the relation between their angular size 
■!? and a number of quasar groups A'' with this size is characterized by a power-law 



2 ! N r^i}-'^ with d w 2.08. 
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The large-scale structure of the Universe is the structure of the galaxy distribution 
in whole observable Universe volume. The regions of higher galaxy number density are 
CN ■ distinguished in the structure. They are filaments of galaxies and galaxy clusters [I] , [2] 

and planes composed of galaxy clusters [5]. These planes are also called superclusters. 
Huge voids are bordered by planes and may be permeated with filaments. Apparently, 
the large-scale structure is sponge- like Scales of these structures equal tens and 
hundreds of megaparsecs (Mpc). Investigation of the large-scale structure of the Universe 
is necessary for understanding of galaxy evolution and physical laws of the Universe 
evolution. 

At present, the geometrical properties of the large-scale structure are investigated 
through galaxies, galaxy clusters [5], [6], [7], [8], [9], quasars |10], [11], [12] and the CMB 
temperature anisotropy [13], [H]. 

The procedure of selection of the main markers of the large-scale structure from the 
general background contains hardly removable mistakes. This is caused by the fact that 
we cannot eliminate mistakes related to projection of foreground galaxies on a cluster, 
overlapping of cluster projections and overlapping of small galaxy groups in a narrow 
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cone of view in a two-dimensional projection of the observable galaxy distribution. Be- 
sides, there is the selection effect when galaxies with apparent magnitude values close to 
some average value are included in a cluster (the cluster is supposed to be composed of 
approximately similar galaxies). 

Investigation of three-dimensional galaxy distribution so far gives not quite reliable 
results due to low quality of sources spectra used for estimation of distances to galaxies 
through the Hubble law (galaxy velocity is proportional to its distance). In addition to 
mistakes relating to peculiar (non-Hubble) galaxy motion, estimation of distance through 
redshift depends on choice of cosmological model in which the distance is computed as a 
function of redshift (redshift-space distortion). 

Therefore, study of statistical and topological properties of the large-scale structure 
through only two-dimensional galaxy distribution in large areas of the celestial sphere 
is relevant. Of course, two-dimensional projection distorts quantitative properties of the 
large-scale structure. However, two-dimensional projection enables to search for the uni- 
versal properties of the galaxy distribution which is typical for the whole observable 
Universe due to high precision of the galaxy equatorial coordinates estimation and large 
areas of modern surveys. The main difficulty of such researches is related to the fact 
that wide-band surveys are composed of narrow-band surveys which differ in precision 
of coordinates estimation and sensitivity of radiation flux measurement. This lead to 
inhomogeneity of resulting survey. Processing of such inhomogeneous survey may lead to 
revealing of fictitious voids and clumps. 

In this paper the processing of available at present data is performed: 

- on the quasar distribution on the celestial sphere (according to the seventh data 
release of the SDSS quasar catalogue [T5]): 

- on the CMB temperature anisotropy [16j. 

The purpose of the analysis is finding of the most general properties of the large-scale 
structure. The main result of the work is revealing of the fractal properties of the quasar 
distribution on the celestial sphere. 

2. Photometrical properties of the SDSS quasars 

The SDSS catalogue is composed of eight surveys of the celestial sphere. In the area of 
the celestial sphere with equatorial coordinates 9^* < a < 16^, 0° < 6 < 55° the catalogue 
contains 105,783 quasars with redshifts 0.0645 < z < 5.4608. 

Areas of these eight surveys overlap. On the one hand, this permits to improve the 
precision of the coordinate estimation in overlapping regions: the survey-averaged errors 
of source's coordinates estimation equal Aa = 0.139", AS = 0.130". On the other hand, 
one can detect more sources in the overlapping regions. Therefore, the number density 
of detected quasars is higher in these areas. So far one cannot find out if these source 
clumps are real or related to inhomogeneity of the SDSS catalogue only. 

We performed the SDSS-quasars distribution analysis in framework of standard cos- 
mological model with parameters: Hq = 70 km ■ ■ Mpc~^ is the Hubble constant, 
Qm = 0.3, = 0.7 are the dimensionless density parameters of dust and A-term re- 
spectively. Cosmological distance to a quasar with redshift z (comoving distance) is 
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determined by the formula: 



c f dz' 

Distances measured in the Hubble distance c/Hq are used below. 

Quasar distribution in right ascension (the equatorial coordinate a) and cosmological 
distance tHq/c is shown on fig. 1. As one can see, there was an epoch of high galaxy 
activity in the Universe evolution. One of the fundamental problems of quasar physics is 
why this epoch was and why there are no quasars in the vicinity of the Local Group of 
galaxies. 




Figure 1: The quasar distribution: cosmological distance - right ascension 



Variation of SDSS-quasars number density in an element of comoving volume of spatial 
cone of the SDSS catalogue is shown on fig. 2 (the number density is normalized to a 
maximal value). The region of higher number density corresponding to the redshift range 
0.35 < z < 2.30 is distinguished. This range is singled out by the catalogue authors after 
taking into account of constant biases in quasar redshift measurements. 

Existence of high galaxy activity epoch indicates that quasars are generally not "stan- 
dard candles" through which the large-scale structure should be studied. However, as- 
tronomers have no other so bright sources. Therefore, we used the five-color photometric 
system of the SDSS catalogue |17j for estimation of SDSS-quasars luminosity function. 

Mean wavelength values of different bands of the system are shown in Table 1. 
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Figure 2: Variation of SDSS-quasars number density in an element of comoving volume 



We have plotted two-dimensional diagrams for different pairs of wavelength bands. It 
has emerged that there is a relatively regular relation only between apparent magnitudes 
in a long- wave part of the photometrical system and Zm- The diagram — Zm is shown 
on fig. 3. 

The regular relation between im and Zm magnitudes holds out a hope that the quasar 
luminosity function obtained for these bands is biased by observation errors in a less 
degree. 

For decrease of significance of evolution effects we have divided the epoch of the redshift 
range 0.35 < z < 2.30 into six layers. The division is performed according to two criteria: 
1) every layer contains comparable quasar numbers; 2) the quasar luminosity function 
plots for different layers (for im and Zm magnitudes) are qualitatively similar. The second 
criterion follows from the following reasoning. The luminosity function describes the 
luminosity distribution of quasars. According to the adopted at present hypothesis, the 
higher quasar luminosity the higher its mass. Therefore, the luminosity function describes 
the mass distribution of quasars. If the mass distributions of quasars are similar in every 

Table 1: Mean wavelength values of different bands of the SDSS photometric system 
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Figure 3: The diagram of i 



Zm relation for SDSS quasars 



layer we may expect that the same common properties of the large-scale structure are 
displayed in every layer. 

The luminosity function equals number density of quasars with absolute magnitudes 
{Mi, Mi + AMj) (for wavelength band i) in a comoving volume of a layer: 



The quasar absolute magnitudes are adduced in the data base [T7|. 

Direct computation by the formula (1) shows that the quasar luminosity function 
plots for chosen layers are qualitatively similar in spite of unavoidable luminosity selection 
(among the distant sources we can see only the brightest ones). 

The SDSS-quasars luminosity functions for and filters are shown on fig. 4 and 5 
respectively. The empirical quasar luminosity function is usually approximated by a two 
power-law function [TH]: 



where a and /3 relate to branches of bright sources faint sources respectively, M* is 
a characteristic magnitude value. The SDSS-catalogue data don't permit to establish 
the luminosity function shape for faint quasars due to selection effects. Therefore, we 
determine the parameter a of the empirical power-law luminosity function for the bright 
sources branch: (M) = 0* ■ l0-0-Ma+^)M ^ 

The parameters of the redshift layers and of the quasar empirical power-law luminosity 
function for im and Zm wavelength filters are given in Table 2. 
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Figure 4: The quasar luminosity function in different redsliift layers for im filter. The mean 
redshift values for every layer are shown in the inner table 
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Figure 5: The same as on fig. 4 for Zm filter 
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Table 2: The parameters of the redshift layers and of the quasar empirical power-law luminosity 
function for im and Zm wavelength filters 



Redshift layer 


(^) 


Layer bounds 


Quasar number 


a {im filter) 


a {zm filter) 


1 


0.55 


0.39 <z< 0.71 


6471 


-2.88 


-2.83 


2 


0.88 


0.71 <z< 1.04 


7145 


-3.07 


-3.11 


3 


1.20 


1.04 <z< 1.36 


9065 


-3.31 


-3.30 


4 


1.52 


1.36 < z < 1.68 


10060 


-3.33 


-3.30 


5 


1.85 


1.68 <z< 2.01 


9420 


-3.32 


-3.24 


6 


2.17 


2.01 < z < 2.33 


4982 


-3.60 


-3.58 



Figures 4 and 5 shows that the luminosity function curves for layers 4, 5 and 6 inter- 
sect. This may be related to either observation errors or quasars physical evolution. For 
example, galaxy mergers occur in regions of higher galaxy number density. This stimulate 
galaxy cores activity and appearance of quasars. Due to this, the mass distribution of 
galaxies changes and the luminosity function shape changes as well. 

3. Angular correlation function and power spectrum of SDSS quasars 

For general description of the quasar distribution we have computed correlation dimen- 
sion, angular two-point correlation function and angular power spectrum for the expansion 
of the observable quasar distribution in spherical functions in each redshift layer. 

The correlation dimension computation methods are given in papers [19], [20], [21], 
for example. 

The correlation dimension characterizes quasar clumping degree and difference of the 
quasar distribution from a homogenous and isotropic one. The dependence of quasar 
number N (r) in a sphere on its radius r for the chosen redshift range 0.39 < z < 2.33 is 
shown on fig. 6. The points distribution on the plot may be described with a power-law 
with a high correlation coefficient value (0.999): 

N{<r)^r'^% (2) 

where the exponent (correlation dimension) equals d^. = 2.17. 

The power-law (2) is usually considered as an indication of fractality of a spatial 
sources distribution. Similar N (r) dependence is typical for galaxies {dc ~ 1.15 -^ 2.25) 
and it is a large-scale structure common law [20], |21] . 
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Figure 6: N - r relation for 0.39 < z < 2.33 



For each layer we have computed the quasar angular correlation function u ("i?) through 
the usually used formula: 



where ^9 is measured angular distance between quasars, is a quasar number in a sample, 
Np is point number in random Poisson distribution in a celestial sphere area of the SDSS 
catalogue, m is a number of quasar pairs with mutual distance ^, nip is a point 
pairs number in random Poisson distribution with the same distance i). angular distance 
on the celestial sphere using in expression (3) is calculated through the formula: 



where 6i, 62, ai and 0:2 are dechnations and right ascensions of two quasars. 

Every layer presents a large-scale structure variant in an epoch corresponding to a 
mean redshift of the layer. If the cosmological principle (statistical homogeneity and 
isotropy of the Universe) is satisfied the correlation functions of different layers must 
be similar irrespective of errors related to overlapping of two-dimensional projections of 
quasar clumps and catalogue inhomogeneity. 

The correlation function shape is revealed to be similar for each of six layers. An 
example of correlation function plot for the layer 1.04 < z < 1.36 is shown on fig. 7. 

We have computed angular power spectrum of SDSS quasars for each layer using 
methods described in [23], [21]. The power spectrum is usually used for comparison of 
large-scale structure formation models with observations because the power spectrum 




(3) 



^ = arccos (cos 61 cos ^2 cos (ai — 02) + sin 61 sin52). 
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Figure 7: The correlation function of SDSS quasars for the layer 1.04 < z < 1.36 



depends on an adopted cosmological model. This fact enables to estimate cosmological 
parameters. 

We describe the quasar distribution on the celestial sphere by a sum of 5-functions: 

1 oo I 

k=l 1=0 m=—l 

The coefficients a]" of the expansion in spherical functions are determined by the 
formula: 

N 

aT^Y^^r* i^k,^k). (5) 

For a catalogue covering only a part of the celestial sphere, estimations 6^ of true 
expansion coefficient values are calculated through the expression: 

(6) 

where the integrals I"^ = J Y/^dQ, jp = J {Yi^'f dQ. For the whole sphere //" = 0, 

— 1. The integrals //" and J/" are computed numerically. 

The angular power spectrum estimation ui is given by the formula: 

9 



where 

mm - — 

m=—l 

is weight-average of coefficients (6) with weights J/". 

An example of SDSS quasars angular power spectrum (7) plot for the layer 1.04 < 2; < 
1.36 is shown on fig. 8. The angular power spectrum shapes for all six layers are similar. 
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Figure 8: Angular power spectrum of SDSS quasars in the layer 1.04 < z < 1.36 



In case of power- law (2), the angular correlation function and the power spectrum are 
power-laws as well: 

CO (1?) ~ 1?^-^, (8) 

ui ~ P--\ (9) 

The least-squares lines plotted on figures 7 and 8 illustrate the power-laws (8) and (9). 

The parameter values of laws (8) and (9) for six layers are given in table 3. The 
first column is layer's number, the second column is the exponent 7, the third one is 
correlation coefficient values for correlation function plot points (fig. 7), the fourth one is 
the exponent 7^, the fifth column is correlation coefficient values for power spectrum plot 
points (fig. 8). 

The layer-averaged exponent values of laws (8) and (9) are close: (7) = 2.08, (7„) = 
1.92. 
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Table 3: Correlation function and power spectrum parameters for six redshift layers 



Layer # 


7 




7n 




1 


2.02 


-0.71 


2.06 


-0.77 


2 


2.23 


-0.85 


2.00 


-0.94 


3 


2.25 


-0.87 


1.84 


-0.90 


4 


2.00 


-0.92 


1.87 


-0.92 


5 


1.86 


-0.94 


1.84 


-0.92 


6 


2.09 


-0.90 


1.91 


-0.92 



The power laws (2), (8) and (9) indicate possibility of fractal properties of the large- 
scale structure |20]. Note that when these power-laws are satisfied exactly in regions 
of high quasar number density (i.e. u {'d) > 1), the equations dc = 3 — ■y and 7 = 7n 
are satisfied [24]. The SDSS-quasars correlation function (3) is less than unity, therefore, 
these equations are not satisfied for obtained values of the correlation dimension dc = 2.17 
and parameters 7 and 7„. 

4. Fractal properties of the CMB temperature anisotropy according to 

WMAP data 

CMB photons show us the Universe as it was at the recombination epoch. The WMAP 
experiment proves convincingly that CMB temperature angular fluctuations exist. 

For analysis of the CMB temperature anisotropy in WMAP experiment the expansion 
of two-dimensional temperature field in spherical functions is used: 

00 / 

1=0 m=-l 

where Oq = Tq is the mean temperature value. 

An integer number of each spherical harmonic fits in the sphere. A multipole 
number / shows how many periods of the harmonic fits in the interval [0,7r]. A positive 
part of a harmonic corresponds to a patch with positive temperature. The patch size is 

of the order of 6* ^ y = — j — . 

The quantity Ci defined through the expression 

{O'lmO^Vm') = Cl5w5mm' 

is the power spectrum. The angular brackets denote averaging over several scannings of 
the celestial sphere. The power spectrum permits to single out the CMB temperature 
anisotropy for different angular scales. 
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The power spectrum averaged over all m modes (their number is 21 + 1) is used in 
investigations: 

1 



m=—l 

Generally, power spectrum values depend on observer's location in the Universe. For 
statistical description of the CMB temperature anisotropy a value averaged over all ob- 
servers (realizations of observations) is desirable. We have one realization only; therefore, 
for the WMAP data processing the "gaussianity hypothesis" is adopted: temperature 
fluctuations are independent and their amplitudes satisfy the Gaussian statistics for all 
observers. In this relative error of the power spectrum is 



Ci 



21 + 1 



The power spectrum dependence on values of mean multipoles is shown on fig. 9. The 
least-squares line is described with a power-law: 



Ci ~ / 



-1.74 



(10) 



Correlation coefficient equals 0.96. This high value implies that the dependence (10) is a 
significant statistical correlation. 



10 



10' 



10 



10% 



10" 



10-% 




10' 



—I ' I ■ ■ ■ ' ' ' 



_J I I ■ ■ ' ■ ' ■ 







10" 



10' 



10 
L 



10 



10 



Figure 9: CMB power spectrum 



The power spectrum of the CMB temperature fiuctuations qualitatively corresponds 
to the power spectrum of the quasar distribution (9): ui ~ Z^^ ^^. 
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Peaks of CMB power spectrum plot are more relief if the function Ci is used. 

271 

What the fractal properties of the large-scale structure (2), (8), (9) and (10) are 
indicative of? The obvious answer is that the properties are consequences of the fractal 
properties of the initial matter density perturbations which further led to star, galaxy and 
cluster formation due to gravitational instability. Now we observe traces of these fractal 
properties through quasars. This interpretation follows from the hypothesis of Gaussian 
(thermal) spectrum of the initial density perturbations in Newtonian approximation [26] . 
However, it is not quite so in the general theory of relativity because Einstein's tensor 
is not invariant under scale transformation of the Riemannian space-time [27], [2S]- If 
the large-scale structure evolution is described by Einstein's gravity theory the fractal 
properties may not conserve, even if the initial fluctuations had the thermal spectrum. 
In paper [26] a cosmological model is suggested in framework of which the observable 
fractal properties may be explained. It has emerged that fractality may follow from the 
fact that matter energies in clumps don't change and are related to each other by a power 
transformation. The simplest consequence of this condition is a case when clumps sizes 
form a geometric progression. In the next section we show that SDSS-quasar clumps 
satisfy this condition. 



5. Wavelet analysis of the quasar distribution on the celestial sphere 

For estimation of quasar clump sizes we used the spherical wavelet transform [29] . 
A spherical wavelet is a function ip & L"^ {S"^) defined on a sphere and satisfying a zero 
mean condition: 



J 1 + cos v 



52 



The spherical wavelet ip is associated with a fiat wavelet ipE (in polar coordinate system) 
through the inverse stereographic projection: 



, ( ^ ^ 
2ipE I 2tan-,v2 



1 + cos § 

If the fiat wavelet is the "mexican hat" 



/27r 

the spherical wavelet equals 



l-2tan2-, 

^{^M = ^^TT T^^^P (-2tan2^). (11) 

V27r l + cosw V 2 1 



The spherical wavelet (11) is an axisymmetric function. 
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Scaling of the spherical wavelet is performed through the formula 

(a^ — 1) cosi/ + (a^ + 1) 

where 

1 ^ 
tan — = - tan — . 
2 a 2 

A displacement is performed through the replacement of the angle '& by the angle 

1?' = arccos (cos ■& cos a + sin i? sin a cos {cp — /3)), 

a and /3 are the parameters of the displacement on the sphere. 

The spherical wavelet transform of function / (i?, 9?) is defined by the formula 

W (a, a, (3)^ J JZJWMf sin ^ d<f, (12) 

where the bar indicates complex conjugation. The integral (12) is composed of positive 
and negative parts. The positive part arises if there arc points under the hat. The negative 
part arises if there are points under the hat brim. Therefore, the wavelet coefficient (12) 
is positive if point number density under the hat is more than number density under the 
hat brim. The coefficient is negative if number density under the brim is more than that 
under the hat. If number density is equal everywhere the wavelet coefficient (12) is zero. 

Thus, wavelet coefficients (12) describe gradients of surface number density of the 
quasar distribution. The wavelet analysis permits to single out quasar clump areas: the 
wavelet coefficients (12) are positive in these areas. Sizes of the areas are of the order of 
the "mexican hat" size a. The more the "mexican hat" size the lower the wavelet analysis 
resolution, because small clump areas and rarefaction areas of the quasar distribution get 
under the hat when size a is large. 

Estimation of true clump sizes is a difficult problem because the wavelet coefficient 
value (12) depends on a clump size, on a number of quasars in the clump and on a quasar 
distribution in the clump. Therefore, there is no definite relation between the size and 
the wavelet coefficient value (12) of the clump. 

For solution of this problem we should ascertain a relation between angular size 'd of 
a clump and size of the "mexican hat" a. The power spectrum Ci (a) for the spherical 
"mexican hat" (11): 

cT (a) = J ilJa,o,oYrdn, 

m=—l 

This power spectrum Ci (a) as a function of multipole moment has one peak at 
value Imax- Therefore, the spherical wavelet (11) singles out clumps with angular ra- 
180° 

dius ~ i^a- One can find out numerically that the angular size related to the size of 

''max 

„ . , „ , , . 180° 180° , 

mexican hat through expression: a — ^ ~ ~ 

y/TT I 



a- 

max 
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Thus, the continuous wavelet transform (12) in which the spherical "mexican hat" (11) 
is used as a analyzing wavelet has resolution of radian for every scaling coefficient a. 

We describe the distribution of N quasars on the celestial sphere through a sum of 
5-functions (4). The wavelet transform of the function (4) equals 



N 



where 



i=l 



exp 



X. 



+ X 



-2-i , (13) 



Xi = tan — 



1 — cos I?,- 



1 — cos i^j cos a — sin i^j sin a cos {ipi — /3) 
\ 1 + cos i^j \ 1 + cos -di cos a + sin -di sin a cos {cpi - (3) 



The map of distribution of the wavelet coefficient (13) values for a = 0.008 radian is 
shown on fig. 10. Quasar clumps (large quasar groups) look like light spots. The more 
quasar number in a clump the higher brightness of a spot (for spots with equal size) . 

W(a.a.P) 
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Figure 10: A fragment of the wavelet coefficient (13) value distribution and an example of a 
selection of a large quasar chain 



For study of fractal properties of large quasar groups we should determine a number 
of groups Nc with a certain angular size 'dc- For this, we need to adopt a criterion for 
determination of group size. As noted above, a wavelet coefficient characterizes a gradient 
of quasar distribution number density. If the wavelet coefficient is negative, there is a 
region of lower quasar number density (rarefaction) in that place. Estimation of quasar 
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group sizes is carried out as follows. We determine positions of local peaks of the wavelet 
coefficients. These peaks usually group in patches of positive wavelet coefficient values: 
one patch may contain from one to some tens of such peaks. A bound of every patch is 
approximated by an ellipse (fig. 10) a major semiaxis of which is considered as a group 
size. For every angular size value -^c one can compute a number of equal large quasar 
groups Nc- 

The dependence Nc {"^c) for the layer 1.04 < z < 1.36 is shown on fig. 11. This is well 
described by a power law (with correlation coefficient value -0,98): 

Nc ~ ^9-'-°^ (14) 

. . . , ■ . ■ ■ ■ . . ■ I T 
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eg 

Figure 11: Nc ~ relation for the layer 1.04 < z < 1.36 

Thus, the large quasar group distribution in angular sizes is the power law (14). This 
distribution indicates that the large-scale is a fractal set and the exponent d ^ 2.08 is its 
fractal dimension [30]. The power law may be a consequence of the fact that large quasar 
groups are self-similar. For example, let their angular sizes form a geometric progression 

^^ = (15) 

where g < 1 is a geometric ratio. Then, a number of quasar groups with sizes 'di>'dn = 
'doq"'~^ satisfies the law (14): 

i=i 1=1 y 
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6. Conclusion 



The main result of the present work is revealing of fractal properties of the large-scale 
structure of the Universe which are described by the laws (2), (8), (9), (10) and (14). 

Note that de Vaucouleurs tried to find a size distribution law of galaxies using the 
first Zwicky's and Abell's catalogues in 1971 |5j. He concluded that galaxy distribution 
differs from a random Poisson one, but he couldn't find a law analogous to (14) due 
to insufficient volumes of the catalogues. The first effective application of the wavelet 
analysis on a plane for revealing of clusters and filaments on a small area of the celestial 
sphere is demonstrated in paper [6j and estimated angular size values of 33 structures 
with angular sizes from 20" to 11° are presented. Using these sizes one can ascertain that 
the law (14) is satisfied for filaments with fractal dimension d 1.16. 

The fractality of the size distribution of quasar groups may be accounted for in frame- 
work of the fractal cosmological model described in paper [26j. Einstein's equations for 
the space-time metric are shown to have a class of solutions in which the metrics are 
related to each other by a simple transformation: a transition from one metric to another 
is just a multiplication by a constant factor (scale transformation or scaling). 

Perhaps, in the epoch of quasars, the large-scale structure was composed of noninter- 
acting spatial regions with metrics related to each other by a scale transformation. Large 
quasar groups discovered in this paper mark these spatial regions. In this case, clumps 
are physically self-similar. Sizes of two groups are related by a scale transformation: 
= Q'^i-i- The whole set of sizes forms a geometric progression of type (15). In this case, 
a number of groups and their size are related by the correlation (14). 

The fractal dimension value d ~ 2.08 for the large quasar group distribution in sizes is 
compared to that of polygonal path of Brownian particle (length distribution of segments). 
This analogy indicates that initial density perturbations from which large quasar groups 
arise, apparently, had a thermal spectrum. 
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